Spin detection at elevated temperatures using a driven double quantum dot 
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We consider a double quantum dot in the Pauli blockade regime interacting with a nearby single 
spin. We show that under microwave irradiation the average electron occupations of the dots exhibit 
resonances that are sensitive to the state of the nearby spin. The system thus acts as a spin meter 
for the nearby spin. We investigate the conditions for a non-demolition read-out of the spin and 
find that the meter works at temperatures comparable to the dot charging energy and sensitivity is 
mainly limited by the intradot spin relaxation. 
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I. INTRODUCTION 

Electron spins in semiconductors and molecular sys- 
tems are good candidates for qubits due to their relatively 
long coherence times. Manipulation of single spins and 
controlled interaction between pairs of spins are essential 
ingredients for quantum information processing. Single 
spin rotation has been demonstrated in electrostatically 
defined quantum dots using the electron spin resonance 
technique.^ Coherent manipulation of a pair of qubits, 
giving rise to entanglement, has also been achieved in 
a semiconductor double dot (DD) device based on fast 
electrical pulses and operating the dots in the spin block- 
ade regimcjS Single spin rotations together with entan- 
glement generation in principle enable universal quantum 
operations. A further operation for a quantum proces- 
sor is spin detection, which is essential for projection of 
the quantum state after computation and read-out of the 
result. This is the main focus of this paper. 

Single spin detection is also important for future spin- 
tronic devices in general, and various electrical and opti- 
cal schemes have been proposed and demonstrated. For 
example Elzerman et al. demonstrated experimentally a 
single-shot read-out of a quantum dot spin using a spin 
to charge conversion technique;^ while Rugar et al. em- 
ployed magnetic resonance force microscopy to probe the 
state of a single spin-- Other schemes for spin read-out 
involve open quantum dots with an inhomogeneous Zee- 
man splitting and closed DD systems which are coupled 
to quantum point contacts. ^'^ It has been shown theoreti- 
cally that the dc-electrical current and shot noise through 
the dots or the point contacts can provide valuable in- 
formation about the spin state, the energy spectrum and 
the relevant decoherence rates. ^'^ However, for optically 
non-active molecular spins no reliable read-out scheme 
exists. 

In this work we consider a single spin (target spin) that 
interacts with the spins of two tunnel-coupled quantum 
dots, as shown schematically in FiglU and demonstrate 
how to probe its state by monitoring the average elec- 
tron occupation on one of the two dots. The interaction 
between the target spin and the spins on the DD induces 



FIG. 1: Schematic illustration of the proposed system for 
a single spin detection. Two tunnel-coupled quantum dots 
are connected to leads enabling current to pass through. A 
nearby spin interacts with the spins on dot 1 and this inter- 
action induces a different Zeeman splitting in the two dots. 
In an external magnetic field and under microwave irradia- 
tion the spins on the two dots can be rotated with a Rabi 
frequency and the average electron occupation exhibits res- 
onances which are sensitive to the state of the nearby spin. 
A charge detector is capacitively coupled to dot 1 and is used 
to monitor a change in the occupation. In the rotating frame 
with respect to the incident microwave field the lead electrons 
have a Fermi energy that depends on their spin state with a 
relative energy difference £o = fii^o where ojq is the frequency 
of the driving magnetic field. Vsd is the bias voltage across 
the dots. The effective Zeeman splitting in each dot, (5i ± J 
and S2 respectively depends on both the orientation of the 
nearby spin and the interaction strength J. 



an effective Zeeman splitting that is different in each dot. 
Also, the sign of the Zeeman splitting depends on the ori- 
entation of the target spin. This target-spin dependent 
asymmetry of the Zeeman splitting makes it possible to 
rotate only one of dot spins and thus results in a target- 
spin dependent lifting of the Pauli blockade. We show 
this lifting of the blockade as a change in the average dot 
occupation, which can be measured by a charge detec- 
tor. Alternatively, the change in occupation is directly 
related to a change in the current through the double dot, 
which, however, might be too small to detect by standard 
dc-measurement techniques. 

In order to realize a non-demolition spin measurement 
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the spins on the DD and the target spin must have dif- 
ferent Zeeman splitting, most hkely to be achieved by 
different ^-factors through g-factor engineering or choice 
of materials. A general spin-spin interaction always con- 
tains so called spin flip terms, which are suppressed when 
the difference in the Zeeman splitting of the dot electrons 
and the target spin is larger than the spin-spin interac- 
tion strength. 

A theoretical investigation has shown that a single 
quantum dot works as a spin meter, but the proposed 
device has the drawback of operating only at low tem- 
peratures, comparable to the energy scale set by the Zee- 
man spittingji In this article we show that the driven DD 
spin detector works at much higher temperatures, com- 
parable to the dot charging energy and we find that the 
main limitation to its sensitivity is intradot spin relax- 
ation. Thus we present a scheme that allows the spin 
state of a nearby spin to be probed noninvasively in a 
single shot, using different electron occupations on the 
double dot for spin up and spin down orientations as a 
basic readout mechanism. 

In the following section we introduce our model and 
discuss technical details of the solution. In Sec. IIIII we 
show how the double dot acts as a spin meter and ex- 
plore how interdot hopping, microwave intensity, temper- 
ature, spin-spin interaction strength and spin relaxation 
influence its performance. We conclude in Sec. II VI bv dis- 
cussing the operation of the spin meter for experimentally 
accessible parameters. 

II. PHYSICAL MODEL 

The total system consists of the DD, a nearby spin, 
metallic leads and a bosonic heat bath. This system is 
modelled by the Hamiltonian 

Htot — Hs + Hieads + Ht + Hb + HsB , 

where Hs models the DD and the spin, HiQ(ifig{^H B^ mod- 
els the leads (heat bath) and Ht{Hsb) models the inter- 
action between the leads (heat bath) and the DD. Specif- 
ically, for the DD and the nearby target spin we write the 
Hamiltonian as 

Hs = Hdd + Hm + Hi, (1) 

where is a Hubbard Hamiltonian describing the 

DD, Hm is due to the applied magnetic fields and Hi 
models the interaction of the nearby spin with the DD 
system. For the DD we have 
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HdD ^^£tn^ - -l^{c\„C2a + cl^Cla) 

.=1 ^ . (2) 

i=l 

that allows up to 2 electrons per dot. The number oper- 
ator is Ui = ^io = Scr cLcicr for dot i = {1, 2} and 



spin cr = {t,i}- The operator c-^ {cia) creates (annihi- 
lates) an electron on dot i with on-site energy e^. 7 is 
the tunnel coupling between the two dots, U is the charg- 
ing energy (intradot Coulomb energy) and V the interdot 
Coulomb energy. The Hamiltonian part due to the ap- 
plied magnetic fields, that breaks the spin degeneracy, 
is 

2 2 
= E T"^' + E cos(a;ot)af , (3) 

i=0 i=l 

where i = Q refers to the target spin and the spin 
operators are defined in the standard way <ri — 
Scrcr' cJcr'''o-cr'CiCT', with (T bclug the vcctor of the 2x2 
Pauli matrices. = gt^sBi is the Zeeman splitting 
due to a static magnetic field Bi along z, and a factor 
gi, fl is the Rabi frequency and luq the frequency of the 
oscillating magnetic field along x. For a single spin the 
oscillating magnetic field rotates the ^-component of the 
spin with frequency fl when A = hujQ. We have ignored 
the effect of the oscillating field on the target spin which 
is a good approximation for narrow-band radiation that 
is only resonant with the spins on the DD (or alterna- 
tively only with the target spin) , a condition that can be 
achieved for example by engineering different g-factors in 
the dots and the target spin. 

Moreover, we assume that the target spin interacts 
only with dot 1, although the basic idea can be extended 
to the most general case when the target spin interacts 
with both dots. As shown below our scheme is still ef- 
ficient provided that the strength of the interaction be- 
tween the target spin and each dot is different, a con- 
dition that is typically satisfied. We consider an Ising 
interaction between dot i — 1 and the target spin i ~ 
of the form 

Hi = ^aoVJ, (4) 

with J being the strength of the interaction. J mainly 
depends on the distance of dot 1 from the target spin 
as well as the actual size of the dot and the target spin. 
Physical values for J for a purely dipolar interaction are 
within the range of a few MHz as shown in Ref. 7. This 
form of interaction is justified when there is negligible 
tunnel-coupling between dot 1 and the nearby spin so 
that to a good approximation hopping can be ignored. In 
addition, spin-flip processes are weak due to the Zeeman 
splitting induced by the static magnetic field and thus 
neglected. This is a good approximation when the differ- 
ence in the Zeeman splittings between the target spin and 
the DD spins is much larger than the interaction strength 
J. Under these conditions the Ising interaction Eq. (U)) is 
a reasonable choice and leads to a non-demolition mea- 
surement. 

The choice of Ising interaction dictates that the com- 
bination of quantum dot system and target electron has 
to be specifically tailored to realise this non-demolition 
measurement. The necessary regime of parameters might 
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FIG. 2: Schematic illustration of the two-electron energy dia- 
gram for the undriven system, luq is the oscillating frequency 
of the driving field. |r+)and |r+) are blocked triplet states, 
\STi) , i = 1, 2, 3 are states with a singlet component. Rl, R2 
and R3 label regions of possible resonances 

be difficult to realise in a gate-defined quantum dot sys- 
tem, e.g. in GaAs, but arises quite naturally in carbon 
nanotube dots probing a molecular spin. For example 
with a typical dipole-dipole spin interaction strength of 
5 MHz, a typical difference in Zeeman splitting of about 
5% and a typical EPR Zeeman splitting of 10 MHz (see 
also Ref. 7) we arrive at a ratio of coupling strength to 
difference in Zeeman splitting of 1/100. Considering the 
flip flop terms as a perturbation to the diagonal Hamil- 
tonian as in Ref. 8, the first order corrections vanish and 
only second order terms contribute, which are typically 
suppressed by a factor of 1/10000. This means that any 
spin flip that can disturb the measurement will take place 
at a much reduced rate, thus making a non-demolition 
measurement possible. 

For the DD with single orbital levels there are in total 
16 states and the maximum number of electrons is 4. The 
many-body states of the DD and target spin system can 
be written in the form \DD)\(tt) where \DD) denotes 
the many-body states of the DD and |(Tt) = | t) or | J, 
) denotes the two possible target spin states. Starting 
with Eq.(IT]) we can define the two uncoupled spin up, 
and spin down, ^ , Hamiltonians in the subspace 
\DD)\ t) and \DD)\ |) respectively depending on the 
state of the target spin. From Eqs.([2]),([3]) and (Hj) it can 
be shown that 

=^DD + cr— H ^ 0-1 + — CT2 

(5) 

-I- ^ hi} cos(cjoi)crf , 

i=l 

with a = -|-1(— 1) for target spin up (down). 

Figure [2] shows schematically the two-electron energy 
level diagram as function of the static magnetic field^ of 
the undriven system (f2 = 0) and for spin up. The upper 
and lower energy levels correspond to the triplets |T-|_) 



and |r_) which as explained in the next section lead to 
spin blockade. The remaining three energy levels, \STi), 
i — 1,2,3, correspond to two-electron states that contain 
a singlet component and thus allow electronic transport. 
For the driven system {ft ^ 0) transitions from the \T±) 
to \STi) states can take place lifting the spin blockade and 
resulting in a change of the average occupation. Within a 
simplified approach if a transition frequency matches the 
frequency of the driving field the corresponding transi- 
tion probability is expected to be high. The regions Rl, 
R2, R3 in Fig. [5] illustrate such a resonant behaviour. 
The magnetic field where the resonances occur depends 
on the energy splitting of the \STi) states and therefore 
on the interdot hopping 7 and the coupling strength J. 
Additionally, as shown in the next section, interference 
effects, for instance when \T±) — IST2), can yield a van- 
ishingly small probability for particular magnetic fields. 

Having described the interaction of the DD with the 
oscillating magnetic field we analyse the dissipative in- 
teractions with the leads and the phonon heat bath. The 
left and right leads are described by a Hamiltonian of 
the form Hieads = J2eka ''ikdlf.^dika, where dj^.^ (dgka) 
creates (destroys) an electron in lead £ — {L, R} with 
momentum k, spin a and energy eik- The interaction 
between the dots and the leads is given by the tunneling 
Hamiltonian 

Ht = '^{tLc\„dLka + tRcl^djika) + H.C., (6) 
ka 

where tL{tji) is the tunnel coupling between dot 1(2) and 
lead L{R) and we consider the symmetric case where = 

To take into account spin relaxation we have consid- 
ered a generic bosonic bath that is modelled as a set of 
harmonic oscillators and is described by the Hamiltonian 
Hb = J2j ^ijoljOij -I- J2j f'^2,jalja2,j- We have as- 
sumed that each quantum dot is coupled to an indepen- 
dent bosonic bath and there are no environment-induced 
correlations between the two dots. The operators a\ ■ 
{cLi.j) create (destroy) a boson in mode j and similarly 
for a\ ■ {0,2. j), while oJij are the corresponding frequen- 
cies of the bath modes. The interaction between the bath 
and the spins of the DD is given by the general model 
Hamiltonian 

HsB = + '^2^^'2-J^2j + H-C, (7) 

j j 

where the spin-flip operators are a~ = cl^Ci^ and 
Ai.j(A2j) is the coupling constant between dot 1(2) and 
the jth mode of the corresponding bath. Hsb allows 
spin-flip processes for electrons in the DD via energy ex- 
change with the bath which, as shown in the next section, 
leads to a leakage current. We consider spin relaxation 
only in the DD since for the target spin an upper limit 
to its relaxation rate is set implicitly by the coupling to 
the leads. In our scheme the relaxation rate of the target 
spin has to be smaller than the electron tunneling rate 
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from the leads to the DD to ensure a measurable change 
in the DD occupations before spin relaxation. 

To investigate the electron occupation of the system we 
employ a master equation approach^ - and derive an equa- 
tion of motion for the reduced density matrix, p, for the 
system of interest that consists of the DD and the target 
spin. The occupation probabilities are given by the diag- 
onal elements of p. It is convenient to eliminate the time 
dependence from the system Hamiltonian Hg and for this 
reason we perform a rotating wave approximation! ^"^i^^ 
This approximation is well-justified only for weak driv- 
ing, i.e., when $7 ^ as in our system. In the rotating 
frame an arbitrary system operator K, is transformed as 
UlKXJz with the unitary operator Uz = cxp{—ia^ujot/2) 

and cr^ = X;Lo<- 

Starting with the total density matrix, Xy s-nd within 

the standard Born and Markov approximations^ we de- 
rive an equation of motion for p by tracing over the leads 
and bosonic bath degrees of freedom, i.e., p=Tte{x} 
where Tr £;{...} means trace over the environmental de- 
grees of freedom. In the rotating frame and having per- 
formed a rotating wave approximation the density matrix 
p satisfies the equation of motion 

p{t) = Csp{t) + CleadsP{t) + CBp{t), (8) 

with the free evolution term 

Csp{t)^-'-[H''s,P{t% 
and the terms due to the electronic leads 

CieadsP{t) = - -^Trsl / drlHrit), 

[U{T)HT{t - T)U\T),p{t) ® Pleads]]}, 

and the bosonic bath 

i^Bpit) = - -^Trel j dT[HsB{t), 

[V(T)HsB{t-T)V\T),p(t)^PB]]}. 

The operators are U{t) — exp[— 1(7^^ -|- Hieads)T/h] and 
V{t) = exp[-i('H5 -I- HB)T/h], with pieads, Pb being 
the equilibrium density matrix for the leads and the 
bosonic bath respectively. The time dependent operators 
are Ci^{t) = Q-f exp(-ia;oi/2), Ciiii) = Ciiexp{+iujQt/2) 
and the Hamiltonian Tig depends on the nearby spin, 
a = 1) for spin up(down), i.e., 

2 

Us = Hdd + y + 2 1 y 2 + 2^ ^cr» : (9) 

i=l 

where we have introduced the magnetic field detuning 
5i = Ai — hijjQ. 

For the numerical calculations we write Eq.([8|) in the 
energy basis. This results in a system of 256 coupled 
equations for all the matrix elements of p which is solved 



numerically taking into account the normalisation con- 
dition for the diagonal elements, X]i=i Pi,i — ^- 
interested in the steady state, pst, that corresponds to 
p = in Eq.®. The quantity of interest is the aver- 
age electron occupation of the DD, for example of dot 1, 
that is calculated as A'^i = Trjnipst}. In the next section 
we present the basic results and explain the influence of 
various system parameters on the average electron occu- 
pation of the DD. 



III. RESULTS AND DISCUSSION 

Before we examine the influence of microwave radiation 
we have to make a choice for the operating regime of the 
DD. DD systems and their physical response are highly 
tunable by adjusting the gate voltages and the source- 
drain bias voltage in the leads. A regime which is easily 
accessible and has attracted a lot of interest is the Pauli 
spin blockade regime which has been demonstrated ex- 
perimentally in various systems such as AlGaAs/GaAs 
and Si/Ge double quantum dot o^^'^'^ as well as carbon 
nanotube dotsJ^ In this regime one electron is confined 
in each dot and the three triplet states are almost equally 
and fully populated. In the absence of spin relaxation 
and microwaves the (1,1) triplet state is blocked from 
moving on to a (0, 2) state by the Pauli exclusion prin- 
ciple {{n,m) denotes a charge state with n{m) electrons 
on dot 1(2)). Thus the electrical current as a function of 
the source-drain bias is suppressed. 

For a fixed source-drain bias in the spin blockade 
regime a change in the occupations of the two dots can 
occur and a microwave-induced current can flow provided 
that the two dots have a different Zeeman splittingJ^ In 
this case the oscillating magnetic field in combination 
with a static field induces coherent spin rotations that 
mix two-electron states and current flows through the 
transport cycle (0,1) ^ (1,1) -J- (0,2) ^ (0,1). When 
the two dots have the same Zeeman splitting the spins in 
the two dots rotate at the same rate in the triplet sub- 
space and therefore the average occupation remains fixed 
and current does not flowJ^ii^ In this case only spin re- 
laxation can give rise to a change in the occupation of 
the dots. 

Inspection of the af terms in Hamiltonian Eq. 12] shows 
that the interaction of the nearby spin with the spins on 
dot 1 induces an effective Zeeman asymmetry between 
the two dots of order J that depends on the orientation of 
the nearby spin. This suggests that a microwave-induced 
change in the occupation of the DD could take place and 
reveal information about the spin state when the dot pa- 
rameters are adjusted to the spin blockade regime. 

A Zeeman asymmetry can in principle arise due to in- 
trinsic factors as in the case where the two coupled dots 
have different g-factors leading to Ai ^ A2 which could 
make the spin detection difficult. In Ref. 17 we have 
shown how to detect a magnetic field gradient and/or a 
difference in the (^-factors in the absence of the nearby 
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FIG. 3: (a) Average occupation of dot 1 as a function of 
the magnetic field detuning ^2 when the target spin is up 
(solid line) and down (dashed line) for the parameters J = 47, 
7 = 1.5 X IQ-^U, ksT = U/WO, hO. = IQ-^U, hwo = IQ-^C/. 
(b) The same as in (a) for J = O.47. 
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FIG. 4: (a) Microwave-induced change in the occupation of 
dot 1 as a function of the interaction strength J for different 
temperatures and when the target spin is up and 7 = 1.5 x 
10"®C/, ha = IQ-^U, hwo = IQ'^U. (b) Maximum difference 
of the spin up and down occupations for the same parameters 
as in (a) with ksT = i7/100. The inset is an enlarged view 
for small J. 



spin that corresponds to J = 0. Within our scheme spin 
detection is efficient when the intrinsic Zeeman asymme- 
try is much smaller than the spin interaction, i.e., when 
|Ai — A2I <C J. Nevertheless, in this work we assume 
for simplicity that Ai = A2 and thus the Zeeman asym- 
metry is due entirely to the presence of the nearby spin, 
and further that J is independent of the applied magnetic 
field. 

For the numerical calculations the internal parame- 
ters of the DD are adjusted to zero energy detuning 



and specifically we choose i?(l,l) — i?(0, 2) = 0, with 
E{n,m) the energy of the charge state {n,m). For a 
practical realization this configuration could be achieved 
via adjusting the gate voltages that define the confining 
potential of the double dot.^^ We choose for the inter- 
dot Coulomb energy V — U/2 and for the on-site en- 
ergies El = ~U/2 and £2 = —U. The bias voltage is 
= (ml — [i'R)l& = U /2e and is applied symmetrically, 
thus /iL — U/A and hr = —U/A, with ^L{^^R) being 
the chemical potential of the left (right) lead. When the 
interdot hopping satisfies 7 ^ C/ and the temperature 
ksT < U/80 the current as a function of source-drain 
bias is suppressed due to spin blockade and each dot con- 
tains a single electron. 

Figure ^a,) shows the average electron occupation of 
dot 1 as a function of the magnetic field detuning 62 for 
the two possible states of the target spin. We consider no 
spin relaxation and therefore we set Aij — A2,j = in 
Eq.([7]). The dot occupation exhibits resonances (peaks) 
due to intradot spin rotations induced by the oscillat- 
ing magnetic field and interdot tunneling, and it is con- 
stant close to unity far from the resonances due to spin 
blockade. For each spin configuration there are in to- 
tal four peaks whose positions depend on 7 and J. The 
two outer peaks correspond to resonances Rl and R3 
(see Fig. ^ and we would expect a third resonance at 
tudQ = (Ai -I- A2 ± J)/^, corresponding to R2, but this 
resonance is split by an antiresonance, resulting in the 
two inner peaks. In terms of spin dependent detuning the 
condition for the antiresonance is |i5i ± J| = \S2\- This 
very symmetric situation together with the microwave 
driving leads to the emergence of an eigenstate of the sys- 
tem with purely (1, 1) triplet components, thus leading 
to spin blockade. It can be shown from the steady-state 
occupations that for this detuning the 1502) = |0,t4'> 
state is unoccupied and the current is suppressed. Away 
from the symmetry point, when k > ^Ml/ J , defining the 
distance from the symmetry point as k = \5i ± J| — |52|, 
the eigenstate is no longer a pure triplet and current can 
flow, resulting in the two inner peaks left and right of 
the antiresonance. We can try to gain an intuitive under- 
standing of the outer peak positions. Naively one would 
expect the outer peaks to appear when either the spin in 
dot 1 is on resonance or the spin in dot 2 is on resonance. 
However, for finite interdot hopping the outer peaks are 
somewhat shifted from the positions that one would ex- 
pect for independent spins since intradot spin rotations 
take place with interdot hopping. As a result the shift 
is large when 7 is large. Interdot hopping leads to delo- 
calisation of the electron spins via the resonant coherent 
transitions | t,!) o |0,t-l> and | 4.,t) 10,14-), with 
an amplitude proportional to 7, that populate the |0, t4-) 
state and thus lead to a change in the populations. In 
addition, the populations of the | t, \.) and | |, t) states 
are unequal leading to a mixing of the (1,1) singlet and 
the 5z = triplet that depends on the magnitude of J 
and 7. When J is small the two outer peaks overlap 
with the inner peaks as shown for example in Fig. [3Ub), 
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FIG. 5: Relative occupation of dot 1 as a function of magnetic 
field detuning S2 and interdot hopping 7 for spin up and down 
for tlie parameters J = 9 x 10"^I7, ksT = U/WO, hO. = 
W~^'U, hujo = W~'^U. The lower frame shows the absolute 
difference of the spin up and down occupations. 



though the important point is that the spin up and down 
signals can still be distinguished (see also below) making 
the measurement feasible. When J — resonances do 
not occur since the spins in the two dots rotate at the 
same rate and the IS'02) state remains unoccupied. 

The change in the dots' occupations is large when the 
Rabi frequency, fl, is smaller than j/h 01 the same order 
of magnitude, a condition that allows interdot hopping 
while intradot spin rotations take place. In the oppo- 
site limit spin rotations dominate and the average dot 
occupation remains close to unity in both dots. If J is 
much smaller than 7 the mixing of two-electron states is 
weak, whereas in the opposite limit interdot transitions 



are suppressed and the microwaves have no significant 
effect resulting in a rather small resonant change. There- 
fore, probing a small J needs a small interdot hopping 
7- 

To quantify this effect we show in Fig.lDJa) the relative 
occupation on resonance as a function of J for a fixed in- 
terdot hopping and for different temperatures when the 
nearby spin is up. A spin down results in the same be- 
haviour. The relative occupation of dot 1, N^, is calcu- 
lated as Nl — Ni — Nf where Nf is the occupation 
on resonance, i.e., the value of the occupation at the 
strongest peak, and is the occupation off resonance 
(background occupation) . The weak feature in the curves 
at J ^ 27 is due to the outer peaks (Rl and R3 regions 
in Fig. [2]) becoming distinguishable from the inner peaks 
(R2 region). For temperatures fc^T < U/80 the back- 
ground occupation is fixed, Nf ^ 1, due to spin blockade 
and hence the relative occupation is essentially tempera- 
ture independent. In this regime there is to good approx- 
imation one electron in each dot in the \T±) states. With 
increasing temperature spin blockade is gradually lifted 
and the background occupation increases. All one- and 
two-electron states acquire a finite population and even 
three-electron states, for instance (I, 2), become occupied 
and have to be included in the dynamics of the density 
matrix. This happens since the lengthening tail of the 
Fermi-Dirac distribution of the lead electrons leads to 
the opening of additional transport channels. The exact 
temperature dependence of depends on various fac- 
tors such as coupling to the leads, spin relaxation rate, 
as well as the applied source-drain bias. Even though 
this dependence may not be monotonic in all cases for 
high enough temperatures {ksT > U) the resonances 
cannot be clearly resolved and ^ 0. From Fig. IH^a) 
we conclude that the DD detector has a higher tempera- 
ture range of operation compared with a single dot, since 
the charging energy is the relevant energy scale. A simi- 
lar increase in operating temperature has been predicted 
for an undriven DD read-out of a charge qubitf^^ In a 
spin read-out situation we are not only interested in the 
height of the resonant peaks, but we want to distinguish 
two target spin states. Thus the figure of merit for a spin 
read-out has to be the maximum difference in population 
for target spin up and target spin down. In Fig. Sfb) 
we plot the maximum difference AiV^" of the spin up and 
down occupations as a function of J and for a fixed inter- 
dot hopping. For J 27 the maximum difference occurs 
at the inner peak of the spin up (down) occupation with 
S2 < {S2 > 0). The results indicate that a large differ- 
ence can be induced making possible the discrimination 
between spin up and down states. 

As shown above the achievable difference in dot occu- 
pation depends on a range of parameters. Figure [S] shows 
a contour plot of the average occupation as a function of 
detuning and interdot hopping for a fixed spin interaction 
strength, 1^ The occupation exhibits a distinct resonant 
pattern for both spin up and down and further it enables 
the two possible outcomes to be distinguished in a range 



7 







T, 


= 2 X lo^^n 


■ ■ ■ X 


= io'"n 


... 


= 2 X lO^^fl 




= 4 X 10"^£2 


Y, 


= 8 X lO^'Sl 



relaxation rate is much smaller than the tunneling rate 
through the double dot and the system weakly deviates 
from the spin-blockade regime. As seen in Fig. [Hlthe ef- 
fect of spin relaxation becomes important as 7s increases 
which in turn leads to a decrease in the relative occupa- 
tion. The sensitivity of the spin detector is limited by 
the minimum detectable change in the occupations. The 
optimum resolution can be achieved when the driving is 
efficient. This happens when the Rabi frequency, which 
is controlled by the intensity of the oscillating field, is 
larger than the spin relaxation rate of the dot spins and 
the target spin, as well as the tunneling rate through the 
DD. 



FIG. 6: Relative occupation of dot 1 as a function of the 
interdot hopping for different spin relaxation rates when the 
target spin is up and the parameters J = 9 x W~^U, ksT — 
17/100, hO, = IQ-^U, hiuo = 10"^(7. A spin down shows the 
same behaviour. 



of interdot hopping. Our calculations confirm that this 
is a robust behaviour that occurs for other values of J 
in the range ^ (10~^ — 10~^)U. However, as explained 
above for ^ ^ J the occupation peaks decrease and this 
could make the spin detection relatively difficult. 

In addition to the temperature effect the background 
average occupation increases due to spin relaxation and 
as a result the microwave-induced resonances cannot be 
clearly resolved since the relative occupation in both dots 
drops. Spin relaxation and decoherence will also infiu- 
ence the peak height of the resonances, since they in- 
hibit coherent spin rotations. Spin-fiip processes which 
take place because of the interaction of the DD spins 
with the bosonic bath described by Eq.© allow incoher- 
ent transitions between two-electron states, for example 
t, t) ^ I i, t), I t, i) which in turn populate the 15*02) 
state, lifting the spin blockade, and thus increasing the 
background occupation. This happens even in the ab- 
sence of the nearby spin, i.e., when J = though spin 
blockade can still be recovered as shown in Ref. 16 de- 
pending on the spin relaxation rate and the coupling to 
the leads. 

To examine the effect of spin relaxation on the driven 
DD spin detector we have calculated the change in pop- 
ulation for various spin relaxations rates 7s . Results 
are shown in Fig. |6] for the relative occupation of the 
strongest (outer) peaks as a function of hopping when 
the nearby spin is up (the same behaviour results for spin 
down). We have taken 7s = ■k\K\^ D{5e)[2n{5e,T) - l]/h, 
with the Bose function n{Se,T) = [exp{6e/kBT) — 1]^^, 
and 5e — hojQ. D is the density of states for the bosonic 
bath that is taken constant and also |A| = |Aij| = |A2j| 
in Eq.Q. This expression for 7s can be derived by as- 
suming a single spin with Zeeman splitting Se coupled to 
a bosonic bath at temperature T which we assume to be 
the same as the temperature in the leads. We focus on the 
most interesting experimental regime in which the spin 



IV. CONCLUSIONS 

In summary, we have suggested an electrical scheme to 
probe a single spin which makes use of two serially tunnel- 
coupled quantum dots connected to metallic leads. The 
spin is located at some distance from the dots, which has 
to be smaller that the typical interdot separation, and 
the total system is in a static magnetic field under the 
application of a microwave magnetic field. The spin in- 
teracts with the spins on the dots and this interaction 
results in an effective Zeeman splitting that is different 
in the two dots. Due to an electron spin resonance effect 
the electron occupations of the dots exhibit resonances 
which reveal information about the state of the nearby 
spin. In particular, the ac-driven DD spin detector pro- 
vides an explicit signal in the induced occupations for 
both spin orientations and enables the spin state to be 
probed noninvasively in a single shot provided that the 
target spin has a different g-factor from the DD system, 
a condition that is typically satisfied. 

We identified a range of parameters for which the sys- 
tem can operate and analysed how we can tune its sen- 
sitivity with the interdot hopping and intensity of the 
microwave field that defines the Rabi frequency for spin 
rotations. The operation of the DD detector depends on 
a lifting of a Pauli spin blockade and therefore it can op- 
erate at much higher temperatures than the single dot 
which is limited to temperatures comparable to the Zee- 
man energyii For instance for a charging energy of 10 
meV and interaction strength in the range ~5 MHz the 
resonances survive up to temperatures of a few tens of 
Kelvin. To achieve a similar operating temperature with 
a single dot, magnetic fields of a few tens of Tesla and cor- 
respondingly microwave frequencies of several hundreds 
of GHz would be necessary, conditions which are available 
only in specialised laboratories. The sensitivity of the de- 
tector is limited by internal spin relaxation which essen- 
tially leads to a small change in the occupations and as 
a consequence the resonances cannot be clearly resolved. 
For an efficient read-out the tunnelling rates from dots 
to leads and the microwave-induced Rabi frequency have 
to be larger than all the relevant spin relaxation rates. 

Finally, the change in population effected by the mi- 
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crowave field has to be seen in relation to the relax- 
ation time of the target spin. On one hand, the spin 
read-out has to be finished within the relaxation time 
of the target spin, otherwise random spin flips will ob- 
scure the result. On the other hand, for a given change 
in population, ANi, of dot i = 1,2 a certain number 
of electrons, n, has to pass through the double dot and 
be counted by the charge detector. To achieve reliable 
statistics the requirement n > 1/ AiV? has to be fulfilled, 
since then fluctuations in the average number are smaller 
than the change in population that we want to distin- 
guish. At a given tunneling rate T through the device 
this determines the minimum time of the measurement 
T™ = n/T = l/{ANfT) < Ti, which must be smaller 
than the spin relaxation time of the target spin. We ar- 
rive therefore at a minimum tunneling rate through the 
dot: r > l/{ANfTi). However, we are not free to in- 
crease the tunneling rate arbitrarily; once the tunneling 
rate approaches the Rabi frequency the resonance peaks 
in the population disappear. By demanding lOF < Q we 
obtain the minimum resolvable change in population as 
AN, - ^/T0/(J\U). 

We conclude by estimating the feasibility of measur- 
ing the state of a molecular spin system, for example 
Sc@C82, which can be coupled with a carbon nanotube 



double quantum dot.''^ Such carbon based systems are 
promising candidates for quantum information process- 
ing and around several Kelvin have a Ti 1 s (Refs. 20, 
and 21). Then with a Rabi frequency of 10 MHz we ar- 
rive at a minimum resolvable change in population of 
AN, = 0.001. An intcrdot hopping of 7 = 10 MHz 
and a spin-spin interaction of J « 5 MHz would lead 
to ANi ~ 0.07 in the case of no spin decoherence. With 
a spin decoherence rate of 7^ = 1 MHz this reduces to 
ANi ~ 0.007, indicating that a single spin read-out with 
realistic parameters at liquid helium temperatures is fea- 
sible. 
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